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ABSTRACT 

We present a comprehensive study of the periodic variations observed in OGLE I-band light 
curves of SMC Be/X-ray binaries, discovering new optical periodicities in 9 systems. We 
find that these periodicities derive from a number of mechanisms, notably disturbance of the 
decretion disk on the orbital period of the system, and aliased non-radial pulsations. We de- 
velop metrics that allow these mechanisms to be distinguished on the basis of the shape of the 
folded optical light curve, and use these metrics to categorise the periodicities present in ~50 
SMC binary systems. We conclude that extreme care must be taken in the interpretation of the 
OGLE light curves since only around 30% of the periodicities present can be unambiguously 
attributed to orbital periods. 

Key words: stars: emission-line, Be - circumstellar matter - stars: oscillations - X-rays: 
binaries - Magellanic Clouds 



1 INTRODUCTION 

Be/X-ray binaries are an important subclass of high-mass X-ray bi- 
nary (HMXB) that consist of a Be star and a neutron star in a binary 
system. Mass transfer is believed to occur initially via a semi-stable 
equatorial decretion disk around the Be star, some fraction of which 
is captured by the neutron star at, or near, its periastron passage. 
These episodes of mass transfer result in X-ray outbursts, during 
which the spin period of the neutron star can often be observed. 
Knowledge of the orbital period of the system is a key diagnos- 
tic tool as it allows the system to be placed on the Corbet diagram 
of Pspin vs P or b ( Corbet 1986) that is used to identify the broad 
classes of HMXB systems. Orbital periods for many Be/X-ray bi- 
naries have been determined from period analysis of the X-ray light 
curves (driven by the spacing of the X-ray outbursts) and Doppler 
variations in their pulse periods, but are still unknown in many sys- 
tems. For this reason, ways have been developed to try to use the 
optical emission as an alternative means to determine the orbital 
period. 

The optical signature of Be/X-ray binaries derives from a 
number of components: (i) the Be star itself, (ii) the Be star decre- 
tion disk, and (iii) any transient accretion disk around the neutron 
star. The optical light curves of Be/X-ray binaries show consider- 
able differences in overall form, but often show long-term large 
variations. Similar, if not identical, effects are seen in isolated Be 
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stars (jMennickent et al. 2002 1, and so these variations are typically 
interpreted as due to changes in the decretion disk size. 

Periodicities, as opposed to aperiodic long-term variations, in 
the optical signature are expected from several sources on several 
timescales: 

(i) Radial and non-radial pulsations of the Be star are often 
present in the 0. 1-2.0 day period range. 

(ii) Perturbation of the decretion disk by the orbiting neutron 
star as it passes periastron, indicating the Be-star orbital period, 
usually in the range 10-500 days. 

(iii) Semi-stable non-uniformities of the decretion disk, such as 
disk warping, present on longer (years) timescales. 

The on-going OGLlQ coverage of the SMC, which effectively 
started with OGLE II in 1996, has provided a wealth of optical (I- 
band) data on the counterparts to the X-ray binary systems. With 
essentially daily coverage for most of each year, it has been possible 
to relate the optical to the X-ray data and frequently explain the 
characteristics of the X-ray behaviour. In addition, regular periodic 
behavioural patterns in optical data have, in several instances, been 
clearly identified with the binary period of the system - see, for 
example, McGo wan et al.| (2008) who showed correlated optical 
and X-ray outbursts every 138d in SXP46.6. However, Schmidtke, 
|Cowley, &" Udalski (2006) have indicated that care must be taken 
when interpreting all modulations as due to binary periods, since it 
is possible for non-radial pulsations in the Be star to beat with the 

1 http://ogle. astrouw . edu .pi/ ogle 3 /xrom/ xrom. html 
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daily sampling regime and produce peaks in power spectra that may 
be confused with binary behaviour. Furthermore, there is evidence 
for possible super-orbital modulation in many of these systems, a 
phenomenon well characterised by the behaviour of the supergiant 
system SMC X-l (see, for example, Trowbridge, Nowak, & Wilms 
(2007 )) that may further complicate the study of these systems. An 
extensive study of such possibilities is presented in Rajoelimanana, 
ICharles, & Udalski|pOTT) . 

In this paper, we present a systematic analysis of the OGLE 
light curves of 49 SMC Be/X-ray binary systems in which we have 
searched for both persistent and transient periodic behaviour. We 
present evidence of both pulsation-related and orbital variations 
visible in the OGLE light curves from a sample of SMC sources, 
often transient in nature, and introduce new metrics that allow some 
differentiation of the various periodic signals in the light curves. 



2 DATA ANALYSIS 

2.1 The OGLE datasets 

The OGLE datasets used in this analysis provide both I-band and 
V-band photometry. The bulk of this work makes use of the I-band 
data because it is both much more comprehensive in sampling, 
and we expect that decretion disk emission contributes a signifi- 
cant fraction of the I-band flux. The V-band light curves have been 
used as supporting data in some analyses. 

Some example OGLE I-band light curves are shown in Fig- 
ure [T| illustrating various examples of optical variability in SMC 
binary systems. SXP0.92 shows an essentially stable I-band flux 
through the 12 years of OGLE II and OGLE III monitoring. 1300 
I-band observations were taken with a mean measurement error of 
5I=9.3mmag; the source varied by only AI=85mmag overall on a 
mean of 1=16.39. SXP2.37 was covered only by OGLE III, and so 
has only -~700 measurement points, but varies by more than one 
magnitude around a mean of 1=14.88. SXP8.80 shows a similar 
variation, in this case covered by both OGLE II and OGLE III mon- 
itoring. SXP756 shows very clear periodic brightening by AI~0.5 
over a stable level of I~ 14.6 every ~350 days. 

Table [TJpresents an overview of the data quality and coverage 
of the I-band and V-band photometry of the 49 objects studied. For 
each light curve, we indicate the date range covered, the number 
of data points, the mean magnitude, the average measurement error 
(S) and the long-term magnitude range (A). Overall, a typical mean 
I-band magnitude is I~15 with a measurement error <5I~7 mmag 
on each observation. 

Although the observations are optimally at daily intervals, ob- 
servation limits mean that a typical year of OGLE monitoring pro- 
duces ~100 usable observations spread over ~220 days, and hence 
the monitoring periods are separated by ~140 day gaps that are 
clearly visible in Figure [T] A typical light curve contains ~350 
points from OGLE II (if available) and ~800 points from OGLE 
III. Figure [2] shows the separations of adjacent data points in the 
light curve of SXP327, which may be considered typical. 

2.2 Analysis method 

We used the fast implementation of the Lomb-Scargle (L-S) pe- 
riodogram (Press & Rybicki 1989) to search for periodic signals 
within the dataset, this method generally being preferred for data 
sets with unequal sampling and data gaps. Periods were searched 
in the range from 2 days (corresponding to the Nyquist frequency 
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Figure 2. Separations of observation times of consecutive observations for 
SXP327, that may be considered typical of the OGLE light curve sampling. 



for nominally daily sampled lightcurves) up to 1000 days, after 
which sensitivity is reduced due to the finite length of the light 
curves. Significance levels for each periodogram were determined 
by using Monte-Carlo simulations wherein the data points were 
randomised in time within the lightcurve, whilst maintaining the 
same overall time structure. A periodogram was derived for 100000 
randomised light curves, and the power level containing 99.99% of 
the peak powers was determined. The errors on the derived periods 
were also derived by Monte-Carlo means, using a bootstrap-with- 
replacement simulation performed with many iterations, the error 
on the period being derived from the distribution of periods recov- 
ered from typically 10000 trials. 

As the Lomb-Scargle method has limited sensitivity to detect 
high-frequency modulation in the presence of low-frequency (i.e. 
long-term) trends in the data, detrending of the light curves was 
performed to remove long-period signals. A rolling-mean calcula- 
tion over a specific time window (typically 50-300d), moved in one 
day steps, was used to determine a detrending model that was then 
subtracted from the light curve. Such detrending acts as a combined 
high-pass frequency filter and mean subtraction that should allow 
the Lomb-Scargle algorithm to work with optimum efficiency when 
searching for small variations on top of larger long-term variations. 

The sensitivity to periodic variations will depend on a number 
of factors, principally the period and amplitude of the signal being 
searched for, the length, sampling and noise levels of the light curve 
being searched, and the ability of the detrending to successfully re- 
move low-frequency noise from the light curve. We performed a 
number of simulations to estimate the sensitivity of our analysis 
to periodic variations. In each case, we created fake light curves 
with a known periodicity but the same sampling, mean flux, and 
random measurement noise as real light curves. These fake light 
curves were then fed through the same analysis process as the true 
light curves to allow us to assess under what conditions the simu- 
lated periodicity is reliably reconstructed. This procedure was car- 
ried out both for 12-year light curves, in order to assess the sensi- 
tivity to a persistent periodic signal, and for 1-year light curves, to 
estimate the sensitivity to transient periodic emission. 

We analysed the basic sensitivity to persistent periodic vari- 
ations at a period of 50d for 4 sample light curves (SXP6.88, 
SXP74.7, SXP214 and SXP327) representing different lengths and 
intrinsic measurement error levels (see Table [TJ. The sensitivity 
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Figure 1. Long-term OGLE II and OGLE III I-band light-curves for SMC sources. The light curves shown many different overall morphologies typical of Be 
systems - see Mennickent et al. 1 2002 1 for further examples. 



limit for full OGLE II + OGLE III light curves was found to vary 
between persistent signal amplitudes of 1.5 and 4mmag depending 
on the light curve quality, and we confirmed that any signal detected 
above the 99.99% confidence limit (at a Lomb-scargle power of 20) 
was reconstructed at the correct period. For one-year datasets, the 
99.99% confidence limit was determined to be at a L-S power level 
of 14, and the sensitivity to periodicities was found to be closely 
related to the intrinsic measurement errors - periodic signals were 
correctly determined for amplitudes above ~5mmag in SXP6.88 
and ~4mmag in SXP214. 

We performed additional simulations to quantify the impact 
of the detrending employed - a series of simulations of a lOmmag 
periodic signal of either 50-d or 100-d period was injected into a 
typical lightcurve, which was then detrended with windows of var- 
ious widths. Only the ratio of the detrending window length to the 
signal period is important, and the loss of L-S power as a function 
of this ratio is shown in Figure [3}. As expected, the L-S power 
is essentially unaffected when the detrending window is signifi- 
cantly longer than the signal period, but the recovered L-S power is 
strongly reduced for periods smaller than the detrending window. 

Finally, to quantify the effect of the sampling pattern and the 
length of the dataset on sensitivity, we performed a set of simula- 
tions in which we created fake light curves matching the sampling 
of a typical composite OGLE-II and OGLE-III light curve, added 
a lOmmag periodic signal of a varying period, and then tested the 
ability of our algorithms to recover the correct period. The sensi- 
tivity was found to be very uniform over a large range of periods, 
with the exception of a small range around one year periods where 
the long-term sampling pattern of the data dominated. While the 
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Figure 3. Recovered L-S power as a function of detrending window:signal 
period ratio. 



correct periods were recovered at all times, the recovered Lomb- 
Scargle power, and hence sensitivity, was reduced by ~40% in nar- 
row windows around one year. This is unlikely to have a significant 
effect for our studies. Unsurprisingly, it also became difficult to re- 
cover the test period when it became comparable with length of the 
light curve. In the case of yearly subsets of the light curves, periods 
beyond ~175 days (around -~75% of the dataset length) could not 
be reliably detected. Using the combined OGLE-II and OGLE-III 
light curve lengths, this corresponds to an upper limit on the pe- 
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Source 


TJD range 


Nj 


I 


.51 


AI 


TJD range 


Ny 


V 


<5V 


AV 


name 








(mag) 


(mmag) 


(mmag) 








(mag) 


(mmag) 


(mmag) 


SXP2.37 


2086.9 


- 4868.6 


707 


14.89 


6.2 


1063 


3326.6 


- 4795.6 


55 


14.80 


4.0 


221 


SXP2.763 


2086.9 


- 4868.6 


71 1 


13.94 


6.0 


192 


3326.6 


- 4795.6 


54 


14.07 


4.0 


57 


SXP3.34 


627.9 - 


4953.9 


837 


15.57 


6.2 


144 


673.8 - 


4954.9 


109 


15.63 


6.0 


129 


SXP6.85 


627.9 - 


4862.6 


646 


14.70 


4.8 


690 


665.9 - 


3712.6 


43 


14.68 


6.8 


322 


SXP6.88 


626.9 - 


4954.9 


1000 


15.37 


6.1 


164 


670.9 - 


4954.9 


99 


15.56 


5.3 


801 


SXP7.78 


2086.9 


- 4952.9 


762 


14.79 


5.1 


196 


3326.6 


- 4953.9 


89 


14.92 


4.0 


75 


SXP7.92 


621.9 - 


4868.6 


1026 


13.82 


4.7 


647 


665.9 - 


4801.5 


90 


13.85 


5.4 


285 


SXP8.80 


466.6 - 


4952.9 


1094 


14.71 


4.8 


952 


466.6 - 


4953.9 


130 


14.71 


4.9 


774 


SXP9.13 


466.0 - 


4954.4 


1077 


16.30 


9.1 


284 


466.6 - 


4954.9 


99 


16.54 


7.5 


126 


SXP11.5 


2203.6 


- 4868.6 


380 


15.26 


6.7 


1002 


3654.7 


- 4801.5 


48 


15.25 


4.0 


596 


SXP15.3 


466.6 - 


4954.9 


1075 


14.47 


4.7 


250 


466.1 - 


4954.4 


93 


14.68 


5.2 


343 


SXP18.3 


466.5 - 


4952.9 


1111 


15.51 


6.3 


1232 


466.6 - 


4953.9 


135 


15.64 


5.3 


862 


SXP22.07 


2104.9 


- 4873.6 


656 


14.05 


6.0 


137 


3315.7 


- 4797.6 


49 


14.21 


4.0 


61 


SXP25.5 


621.8 - 


4954.9 


1067 


15.72 


7.1 


372 


645.9 - 


4954.9 


94 


15.83 


6.3 


400 


SXP3 1 .0 


2104.9 


- 4873.6 


656 


15.31 


6.0 


197 


3315.7 


- 4797.6 


49 


15.49 


4.0 


28 


SXP34.08 


2086.4 


- 4954.4 


724 


16.94 


13.3 


105 


3326.6 


- 4954.9 


69 


16.84 


7.6 


48 


SXP46.6 


2086.4 


- 4954.4 


1857 


14.74 


5.7 


51 1 


3326.6 


- 4954.9 


148 


14.80 


4.0 


310 


SXP59.0 


626.9 - 


4954.9 


1560 


15.41 


6.4 


657 


670.9 - 


4954.9 


140 


15.44 


5.0 


365 


SXP65.8 


621.9 - 


4862.6 


589 


15.73 


6.8 


320 


673.8 - 


3379.6 


39 


15.77 


5.9 


118 


SXP74.7 


466.5 - 


4954.9 


1715 


16.68 


11.2 


179 


467.1 - 


4954.4 


176 


16.89 


7.8 


167 


SXP82.4 


466.6 - 


4952.9 


1105 


15.08 


5.2 


845 


466.6 - 


4953.9 


131 


15.21 


4.8 


577 


SXP9I.1 


2086.9 


- 4865.6 


713 


14.78 


5.3 


126 


3326.6 


- 4791.6 


50 


15.02 


4.0 


86 


SXP101 


2086.9 


- 4954.9 


737 


15.61 


6.9 


127 


3326.6 


- 4954.9 


71 


15.70 


4.3 


67 


SXP138 


2086.9 


- 4952.9 


762 


16.11 


7.8 


294 


3326.6 


- 4953.9 


89 


16.18 


5.2 


125 


SXP140 


2086.9 


- 4954.9 


724 


15.52 


6.6 


550 


3326.6 


- 4954.9 


67 


15.76 


4.6 


371 


SXP152.1 


2086.4 


- 4954.4 


718 


15.43 


6.3 


182 


3326.6 


- 4954.9 


67 


15.69 


5.4 


125 


SXP169.3 


2086.4 


- 4865.1 


712 


15.38 


6.0 


192 


3326.6 


- 4791.6 


50 


15.57 


4.7 


96 


SXP172 


466.6 - 


4954.9 


1075 


14.42 


4.7 


484 


466.1 - 


4954.4 


95 


14.47 


5.1 


169 


SXP202A 


621.9 - 


4954.9 


1040 


14.65 


5.4 


351 


665.9 - 


4954.9 


105 


14.75 


4.8 


184 


SXP202B 


621.4 - 


4868.1 


911 


15.34 


6.0 


259 


670.4 - 


4795.1 


88 


15.52 


5.5 


181 


SXP214 


466.5 - 


4865.6 


1060 


15.22 


4.2 


447 


673.8 - 


4954.9 


109 


15.63 


6.0 


129 


SXP264 


621.8 - 


4954.9 


1067 


15.90 


7.6 


283 


645.9 - 


4954.9 


94 


16.09 


6.5 


146 


SXP280.4 


2086.9 


- 4954.9 


728 


15.38 


6.1 


244 


3326.6 


- 4954.9 


69 


15.56 


4.8 


198 


SXP293 


621.9- 


4868.6 


1037 


14.57 


4.7 


288 


665.9 - 


4795.6 


93 


14.87 


5.1 


143 


SXP304 


627.4 - 


4954.4 


1052 


15.58 


6.1 


68 


665.9 - 


4954.9 


102 


15.75 


5.9 


55 


SXP323 


621.3 - 


4954.4 


1060 


15.24 


5.6 


216 


466.6 - 


4954.9 


97 


15.42 


6.0 


261 


SXP327 


2086.9 


- 4952.9 


762 


16.75 


11.6 


442 


3326.6 


- 4953.9 


89 


16.70 


5.4 


213 


SXP342 


2086.4 


- 4954.4 


1116 


15.12 


6.4 


967 


3326.6 


- 4954.9 


1 16 


15.16 


4.1 


570 


SXP348 


627.9 - 


4953.9 


688 


15.10 


5.2 


963 


665.9 - 


4954.9 


62 


15.08 


5.9 


629 


SXP455A 


627.9 - 


4954.9 


1762 


15.68 


6.5 


718 


665.4 - 


4954.4 


193 


15.63 


5.3 


338 


SXP504 


621.3 - 


4868.1 


1014 


14.79 


5.2 


83 


670.9 - 


4795.6 


84 


15.01 


5.1 


69 


SXP565 


628.4 - 


4948.4 


845 


15.83 


7.0 


224 


665.9 - 


4948.9 


78 


16.02 


7.0 


110 


SXP645 


621.8- 


4868.6 


1030 


14.66 


5.7 


916 


670.9 - 


4795.6 


87 


14.69 


4.8 


132 


SXP701 


621.8- 


4868.6 


1030 


15.70 


7.1 


148 


670.9 - 


4795.6 


87 


16.02 


6.8 


139 


SXP726 


627.4 - 


4953.4 


1022 


15.49 


5.8 


128 


673.8 - 


4954.9 


127 


15.66 


5.7 


135 


SXP756 


466.6 - 


4954.9 


1080 


14.64 


4.7 


718 


466.1 - 


4954.4 


98 


14.91 


5.1 


402 


SXP893 


466.5 - 


4954.9 


1080 


15.95 


7.7 


128 


466.1 - 


4954.4 


98 


16.27 


6.7 


177 


SXP967 


2104.9 


- 4868.6 


690 


14.30 


6.0 


633 


3326.6 


-4801.6 


57 


14.35 


4.0 


147 


SXP1323 


627.4 - 


4953.4 


1021 


14.57 


5.7 


153 


673.8 - 


4954.9 


127 


14.65 


4.7 


124 



Table 1. Basic parameters of 49 OGLE light curves for SMC Be/X-ray binaries. N/ is the number of data points in the I-band light curve. I is the mean I-band 
magnitude. 51 is the mean measurement error on the data points for each light curve. AI is the overall range of brightness in the light curve (in magnitudes). 
TJD = JD - 2440000.5. Equivalent parameters are given for the V-band light curves. 



riods we can detect of ~2000 days, but in practice the detrending 
windows used made this limit much lower. 



3 PERIOD SEARCHING 

The detrended light curves were searched for periodic signals with 
periods between 2 days and 1000 days looking at both the whole 
light curve and the annual campaigns, subject to the sensitivity lim- 



its discussed in the previous section. In the following sections, we 
report frequencies found in the formally acceptable region from 2 
-1000 days, and reserve any comment on the true origin of these 
periodicities until later sections. Table [^indicates the periodic sig- 
nals detected. We note also the periods previously published in both 
optical and X-ray light-curves for comparison. 

We emphasise here that the recovered period in the 2-1000 
day region may not always represent the true period, but may result 
from aliasing of a periodic signal from outside the formal search 
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Source 


X-ray P orb 


Pub P opt 


Our P opt 


Detrend 




(days) 


(days) 


(days) 


(days) 


SXP 2.37 


18.38l 16 l 


18.62 I 18 ' 


9.30±0.01 


51 


SXP 2.76 




82. lPI 


81.81±0.06 


101 


SXP 3.34 






10.7248±0.0006 


51 


SXP 6.85 


2i 9(16] 


H4[9] 






SXP 7.78 


44.92M 


44.61 3 ' 


44.93±0.01 


101 


SXP 7.92 




36.8I 10 ] 


35.61±0.05 


101 


SXP 8.80 


28.5I 1 * 16 ] 








SXP 9.13 




40.17W 






SXP 11.5 


36.3' 12 ' 








SXP 15.3 




75. lW 


74.32±0.03 


101 


SXP 18.3 


17.7W 


17.7(13] 


17.79±0.03 


101 


SXP 22.1 






83.7±0.1 


101 


SXP 25.5 




22.50' 20 ' 


22.53±0.01 


51 


SXP 31.0 




. fol 

90.4l 2 l 


90.53±0.07 


101 


SXP 46.6 


137.36W 


137P1 


137.4±0.2 


165 


SXP 59.0 


122.25W 


60.21 s ) 






SXP 65.8 




no! 14 ! 


111.04±0.15 


201 


SXP 74.7 


33.38l 16 l 


33.41 4 - 5 ] 


33.387±0.006 


101 


SXP 91.1 




88.21* 5 ' 


88.37±0.03 


101 


SXP 101 




21.9P] 


21.949±0.003 


101 


SXP 140 




19718] 






SXP 169 


68.5W 


67.6PI 


68.37±0.07 


101 


SXP 172 


70.4' 17 ' 


69. 9 [8] 
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101 
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301 
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301 


SXP 701 
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SXP 756 




11.4119] 
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389.9W 


39412] 


393.6±1.2 


501 


SXP 893 






3.7434±0.0005 


51 


SXP 967 




101.94P 1 ] 


101.4±0.2 


151 


SXP 1323 






7.9101±0.0004 


25 






26.ll 8 ] 
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Table 2. Sources for which significant periods have been detected. All units 
are days. References for the published orbital periods are: (1) Galac he et al.| 
(2008); (2)|Schmidtke, Cowley, & UdaTsT^{2006); (3) |Coe et 31.^2005); 
(4)|Edge|(2005); (5) |Schmidtke & CowleylpOOS); (6) |Schmidtke et al] 
(2004); (7)|McGowan et al.N2007^; (8)|Schmidtke & CowleyH2006); (9 ) 
|McGowan et al.p008> ; (10) |Coe et aE]p009> ; (12) |Townsend et al.H2009) ; 
(13)|Schurch et al.|j2009> ; (14)|Schmidtke & Cowley] (2007) ; (15)|Coe et| 



al. 



1 2008 1; (16) Townsend et al.H2011) ; (17) |Schurch et aI.H20H) ; (18) 



Schurch, Udalski, & Coe (200 8); (1 9) |Coe & Ed ge 1 2004 1; (20) Rajoeli- 



manan a, Charles, & Uda lski (20111; (21) Schmidtke, Cowl ey, & Udalski 
20091 



range with other periodicities present in the dataset - usually the 
sampling frequencies of the data. Of particular relevance in the con- 
text of this work is the effect of periods around 1 day, that when 
combined with the sampling pattern associated with the OGLE 
project datasets (Figure [2| can produce complex aliasing patterns 
(often referred to as window functions). Figure [4] shows a L-S pe- 
riodogram of a simulated sine wave with a period of 0.8 days corn- 




Frequency (c/d) 

Figure 4. Lomb-scargle periodogram for 'fake' sinusoid at period of 0.8 
days (/=1.25c/d) superimposed on the SXP1323 light curve sampling pat- 
tern. The sampling frequency (f s ) and Nyquist frequency (/jv) are indi- 
cated by vertical lines. 



bined with a typical OGLE sampling pattern. Apart from the fun- 
damental frequency at 1.25 c/d, a complex set of aliases spanning 
from to ~8 cycles/day, together with the usual reflections about 
the Nyquist frequencies can all be observed. The broad envelope 
of the peaks matches, and is clearly driven by, the data sampling 
shown in Figure|2] 

Although the formal Nyquist frequency of the datasets with 
nominal daily sampling is 0.5 cycles/day (a 2 day period) we can- 
not ignore the effect of periods below 2 days when searching. Due 
to aliasing, periods of very near to 1 day can inject power into the 
periodogram far from the Nyquist frequency, and into the typical 
range for Be/X-ray orbital period searching. For example, a period 
of 0.98d corresponds to 1.02 c/d, creating an alias at 0.02 c/d that 
could easily be misinterpreted as a 50d period. This problem is il- 
lustrated in Figure [5] which shows the detected period (in the for- 
mally useful region between 2 and 100 days) resulting from super- 
position of periodic signals of various periods onto a typical OGLE- 
II and OGLE-III sampling pattern. There is a clear degeneracy in 
the solution as a recovered period could result from a number of 
true periods, and some a priori knowledge of the expected periods, 
or some other discriminator is required to break this degeneracy. A 
possible discriminator, the shape of the folded light curve, is dis- 
cussed in the following section. One further implication of this pat- 
tern is for drifting periods (see Section |53) - while the alias of a 
drifting period will still manifest itself as a drifting period of >2 
days, neither the period or rate of change of the original drifting 
signal can be uniquely determined as a result of this degeneracy. 



4 FOLDED LIGHT CURVE PROFILES 

Where a clear periodicity was detected in a light curve, or annual 
section of a light curve, we created a folded light curve for inspec- 
tion. It was noted that the folded light curves fell into two broad 
categories and were either sinusoidal variations, or sharply peaked 
variations often with 'fast rise exponential decay (FRED)' shaped 
profiles. 

In order to quantify the shapes of the folded light curves, we 
defined two metrics for the folded light curves. The first is the full- 
width half-maximum in phase units, hereafter the phase span (PS). 
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Simulated (input) period (d) 

Figure 5. Recovered period in the 2 - 1000 day region for simulated peri- 
odicities in the range 0.4 - 100 days. 



For a sinusoidal variation, the expected value for this metric is 0.5; 
values below 0.5 would be expected for sharper peaks. The second 
measures the asymmetry of the peak and is defined as the ratio 
between the distances in phase space between the peak and points 
at 10% of maximum on the right and left sides of the peak, hereafter 
the phase asymmetry (PA). For a sinusoid, with a symmetric peak, 
the expected value is 1 ; values greater than one indicate a peak with 
a faster rise than fall. 

The folded light curve analysis was performed for those 
sources exhibiting strong persistent or quasi-persistent (seen in 
more than one year of data) periodicities. The resulting PS and PA 
metrics are listed in Table[3] together with information on the sub- 
sets of data they have been extracted from. 

Two folded light curves are shown to illustrate the principle, 
representing rather extreme cases. The light curve for SXP756 (Fig- 
ure |6j when folded on the orbital period shows very sharp peaks 
occupying a narrow phase range, and showing a large asymmetry. 
The darker shaded region indicates the measured phase span (here 
PS=0.1) and the lighter shaded region indicates the region used in 
calculation of the asymmetry (here PA=6.33 indicates a very asym- 
metric peak). 

The folded light curve for SXP1323 folded on one of the two 
detected periods (Figure[7]l shows much more sinusoidal behaviour 
with PS=0.43 and PA=1.08. 

The folded light curve parameters for those SMC binaries for 
which folded light curve analysis was possible are shown in Fig- 
ure [8] It is apparent that while some of the SMC systems are clus- 
tered around (PS,PA)=(0.5,1.0) point where sinusoids are expected 
to lie, others show a broad range of folded light curve parameters. 



5 DISCUSSION 

The "Match" column in Table ^indicates those sources for which 
the period detected in optical monitoring is also seen in long-term 
X-ray light curves. In Be/X-ray systems, accretion and the result- 
ing X-ray emission normally only occurs when the neutron star is 
close to the Be star decretion disk, and there is no strong evidence 
for stable accretion disks forming around the neutron star. Hence, 
since the X-ray period can only realistically be due to orbital mo- 
tion, we can use this to select a subset of sources where we believe 
the optical period is indicative of the orbital motion. These sources 




Figure 6. Folded light curve analysis for SXP 756 (folded on the 393.6d 
period). The dashed horizontal lines indicate 10% and 50% of the maximum 
value, where the phase span (darker shaded region) and phase asymmetry 
(lighter shaded region) are evaluated. 



1.08 




1.0 
Phase 



Figure 7. Folded light curve analysis for SXP1323 (folded on the 26.17d 
period). The dashed horizontal lines indicate 10% and 50% of the maximum 
value, where the phase span (darker shaded region) and phase asymmetry 
(lighter shaded region) are evaluated. 



are SXP18.3, SXP46.4, SXP169.3, SXP172, SXP504, SXP565 and 
SXP756. There is a strong possibility that SXP280.4 also falls into 
this category, as the X-ray and optically derived periods are a factor 
of 2 different. Where the folded light curve can be parameterised, 
the majority of these sources fall far from the (0.5,1.0) point in Fig- 
ure[8]and so are clearly associated with FRED-like profiles. 

Conversely, those sources where there is evidence that the 
optical period is of non-orbital origin cluster around the (0.5,1.0) 
point. For example, SXP1323 has two strong periods, which is dif- 
ficult to explain as orbital motion. Furthermore, the two observed 
periods are strongly incompatible with being orbital based on the 
Corbet P rb-Ps P in relation for a source with P sp i n =1323s ( |Corbet| 
1986). Similarly, the 11. 4d period in SXP756 cannot be orbital, as 
the 393d orbital period is established beyond doubt. 

We therefore conclude that this 'FRED-iness' analysis and di- 
agram is a useful tool to identify optical periods that are likely to 
be associated with orbital perturbations of a decretion disc based 
on the shape of the folded optical light curve. Both phase span and 
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Source 


Date range 
(TJD) 


Period 


Match 


DW 
fdl 

KM J 


PS 


PA 


Amplitude 
(x 10 -5 Jy) 


Amplitude 

( mmap'1 
\ iiiiiia &/ 


SXP 2.37 


4191-4549 


18.60 


Y 


51.0 


0.24 


3.00 


16.0 


37.9 


SXP 2.763 


all 


81.8 


7 


101 


0.35 


1.38 


14.1 


16.9 


SXP 3.34 


all 


10.72 


7 


51 


0.52 


0.85 


7.5 


38.2 


SXP 7.78 


all 


44.93 


Y 


101 


0.3 


1.5 


14.5 


38.4 


SXP 7.92 


all 


35.61 


? 


101 


0.63 


1.33 


24.0 


23.5 


SXP 15.3 


all 


74.34 


7 


101 


0.53 


0.75 


7.2 


15.1 


SXP 18.3 


3101 -3824 


17.79 


Y 


51 


0.33 


0.86 


21.8 


116.0 


SXP 22.1 


3101 -5000 


84.9 


7 


121 


0.47 


0.88 


11.0 


15.0 


SXP 25.5 


all 


22.53 


7 


51 


0.27 


1.50 


3.6 


21.6 


SXP 31.0 


all 


90.53 


7 


101 


0.33 


5.5 


8.8 


39.4 


SXP 46.6 


all excl. 3100-3800 


137.38 


Y 


165 


0.13 


2.50 


10.9 


25.7 


SXP 65.8 


all 


111.04 


7 


201 


0.33 


2.67 


5.1 


29.4 


SXP 74.7 


all 


33.38 


Y 


101 


0.33 


1.0 


1.9 


30.7 


SXP 91.1 


all 


88.37 


7 


101 


0.33 


1.1 


20.0 


58.6 


SXP 101 


all 


21.95 


? 


101 


0.40 


1.29 


5.9 


31.9 


SXP 169.3 


all 


68.37 


Y 


101 


0.24 


4.0 


7.2 


34.4 


SXP 172 


all 


68.785 


Y 


81 


0.46 


2.33 


8.1 


14.0 


SXP 202B 


all 


224.57 


Y 


501 


0.23 


2.7 


6.3 


28.8 


SXP 214 


all 


4.58 


7 


101 


0.5 


0.6 


13.0 


67.0 


SXP 264 


all 


49.11 


? 


101 


0.2 


2.5 


4.7 


36.3 


SXP 280.4 


all 


127.6 


7 


201 


0.37 


2.33 


7.0 


33.4 


SXP 293 


all 


59.72 


7 


101 


0.37 


1.63 


41.4 


108.0 


SXP 327 


all 


45.9 


? 


101 


0.23 


2.4 


14.0 


205.8 


SXP 455 


all 


74.55 


? 


101 


0.3 


0.88 


- 


- 




all 


209.7 


? 


301 


0.65 


1.11 






SXP 504 


all 


270.1 


Y 


501 


0.4 


3.6 


4.4 


12.5 


SXP 565 


all 


152.4 


Y 


301 


0.2 


1.29 


3.4 


24.6 


SXP 756 


all 


393.6 


Y 


501 


0.1 


6.33 


80.0 


228.5 




non-outburst 


11.404 


N 


31 


0.57 


1.37 


9.7 


25.3 


SXP 893 


1283-3101 


3.74 


? 


51.0 


0.59 


1.0 


2.9 


26.5 


SXP 967 


2200 - 4700 


101.37 


? 


151 


0.33 


2.17 


22.0 


34.4 


SXP 1323 


all 


7.91 


? 


25 


0.55 


0.67 


16.0 


32.9 




all 


26.17 


7 


101 


0.43 


1.08 


24.0 


49.8 



Table 3. Folded light curve parameters for periodic sources in the SMC sample. The match parameter indicates when the X-ray and optical periods match thus 
confirming an orbital period, while a ? indicates that the X-ray period is unknown. Columns indicate the detrending window (DW), the phase span (PS) and 
the phase asymmetry (PA) from folded light curve analysis - see text for more details. 



asymmetry of the outbursts from disk disruption may be, to some 
extent, dependent on several orbital parameters - notably period 
and eccentricity. If we make the naive assumption that the decre- 
tion disks in Be/X-ray binary systems all recover with a similar 
timescale, the phase span parameter will only be effective when the 
disc recovery time is significantly less than the orbital period. Intu- 
itively, more eccentric systems will show a wider range of interac- 
tions between neutron star and decretion disk, while a completely 
circular system would give a much reduced optical modulation (due 
only to radial velocity shifts of the spectrum through a station- 
ary filter). Overall, the method should be most sensitive for long- 
period, eccentric systems, but even reasonably low-eccentricity 
systems should provide a recognizable orbital signature in the opti- 
cal light curve. In summary, then, a FRED-like folded optical light 
curve profile is a strong indicator that the period is orbital. On the 
other hand, a very sinusoidal folded light curve profile is often as- 
sociated with non-orbital periodicities, but could still be indicative 
of orbital motion in a very low-eccentricity system. 



5.1 New orbit determinations 

We have detected previously undiscovered periodicities in the 
OGLE II and III light curves for several sources, and confirmed pre- 



vious optical modulations in many others. In some of these cases, 
the folded light curve analysis gives a strong indication that the 
periodicity is derived from the orbit of the system, allowing us to 
propose them as orbital period determinations. 

Specifically, the folded light curve analyses for SXP31.0, 
SXP65.8, SXP264 and SXP327 provide strong evidence that the 
measured periods at 90.53, 111.04, 49.11 and 45.9 days respec- 
tively, are indeed the orbital periods of these systems. Furthermore, 
the folded light curves of SXP280.4 and SXP967 show quite con- 
vincing evidence of non-sinusoidal (i.e. orbital) modulation at pe- 
riods of 127.6 and 101.4 days respectively (see Figures|9]and 10 1. 
The latter is a confirmation of the detection in Schmidtke, Cowley, 
|& Uda lski (2009), but we note that this is a complex lightcurve, 
difficult to detrend, and the orbital outbursts are only visible in a 
subset of the data. 

SXP 2.37 (SMC X-2) is a particularly interesting case, which 
has been the cause of considerable speculation as to the origin of 
the periodicities seen in the light curve. Previous studies ( |Schurch,1 
Udalski, & Coe 2008 1 have seen periodicities at ~ 6.2, 9.3 and 18.6 
days, which could be interpreted as an orbital period of 18.62d 
with harmonics. However, this interpretation was challenged by 
Schmidtke, Cowley, & Udalski (2009) who interpreted the signals 
as two aliased pulsations at 0.8592d and 0.9008d with the 18.62d 
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Figure 8. Folded light curve metrics for the SMC Be/X-ray binary sample. The solid lines indicate the approximate boundary between regions occupied by 
sinusoidal (likely aliased pulsations) and FRED-like (likely orbital) light curve profiles. 



period deriving from the beat frequency between these two pulsa- 
tions. However, our analysis (Figure [TT) shows that the transient 
18. 6d period is firmly in the 'orbital' region of Figure[8] and we 
thus conclude that this is by far the simplest explanation for the 
three periodicities seen in the periodogram. It has also now been 
shown independently from the X-ray data (Townsend et al.|2011) 
that fitting of the orbital Doppler shifts of the pulsar timing infor- 
mation during a long outburst yields an estimate of the orbital pe- 
riod of 18.38±0.02d. 

For the first time, we confirm that the 152.4d orbital period 
of SXP565, determined from X-ray measurements (Gal ache et al.| 
|2008[ >, can also be seen in the optical light curve (Figure [TTj>, 

5.2 Sinusoidal variations 

The periodicities observed in SXP3.34, SXP7.92, SXP15.3, 
SXP455 (209.7d period), SXP893 and SXP1323 do not show the 
characteristics of orbital modulation of a decretion disk, but instead 
have folded light curves that are essentially sinusoidal in shape. 
Those in SXP3.34 and SXP893 (Figure p~2]> are newly discovered 
periodicities, but for obvious reasons we do not propose them as 



new orbital periods. In particular, the period for SXP893, at 3.74 
days, is incompatible with the usual placement of Be/X-ray bina- 
ries in the Y or b-P ap i n plane. 

SXP1323 is an informative case as it shows two strong and 
persistent pulsations which manifest themselves as 7.91 and 26.17d 
periods in the OGLE light curves. These periods have been at- 
tributed ( Schmidtke & Cowley 2006 ) to the aliasing of non-radial 
pulsations at periods of 0.96 and 0.88 days in the Be star with the 
Id sampling of the dataset. Thus apparent periods at ~8 and ~26 
days, with sinusoidal folded profiles, result from non-radial pulsa- 
tions of the Be star in this system, and this is a plausible explanation 
for at least some of the other non-orbital periods determined. 

As a further example, the second period in SXP756, of 11. 4d, 
is only visible outside of the sharp orbit-induced increases in 
brightness accompanying the Type I X-ray outbursts that are clearly 
separated by the 393. 6d orbital period. This secondary period has 
also been attributed to pulsations of the Be star or an orbiting distur- 
bance in the outer decretion disk, and furthermore is seen to change 
slightly with time ( |Cowley & Schmidtke|2003j|Coe et al.]2003) 
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Figure 9. Periodiograms (left) and folded light curve analyses (right) for SXP31, SXP65.8 and SXP264, confirming previous optical periods, and strongly 
suggesting orbital origins for the detected periods. Dashed horizontal lines on periodograms indicate 99 and 99.9% confidence limits. Dashed horizontal lines 
on folded light curves indicate 10% and 50% of the maximum value, where PS and PA are evaluated (see text for more details). 
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Figure 10. Periodiograms and folded light curve analyses for SXP280.4, SXP327 and SXP967, confirming previous optical periods, and strongly suggesting 
orbital origins for the detected periods. Dashed horizontal lines on periodograms indicate 99 and 99.9% confidence limits. Dashed horizontal lines on folded 
light curves indicate 10% and 50% of the maximum value, where PS and PA are evaluated (see text for more details). 
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Figure 11. Periodiograms and folded light curve analyses for SXP2.37 and SXP565. The first indicates the previously optical period at 18.6 days (now 
confirmed in the X-ray band) is of orbital origin, the second is the first optical detection of the previously detected X-ray period at 152.4 days. Dashed 
horizontal lines on periodograms indicate 99 and 99.9% confidence limits. Dashed horizontal lines on folded light curves indicate 10% and 50% of the 
maximum value, where PS and PA are evaluated (see text for more details). 



5.3 SXP264 - a Rosetta Stone? 

SXP264 is a strange but revealing object because, like SXP756, the 
optical light curve exhibits periods in addition to a well-established 
orbital modulation. Previous studies (Edge 2005) have indicated 
that it shows a ~49d orbital period in optical data. Our time- 
resolved analysis, however, shows a slightly different and more 
complex story. While the expected orbital period is stably detected 
for the first 10 years of OGLE-II and OGLE-III monitoring, the 
peak in the time-resolved periodogram shifts significantly in the 
final 2 years of OGLE-III measurements. 

A deeper analysis is summarised in Figure [13] The first 10 
years of data yield a folded light curve that shows a clear FRED- 
like orbital profile at the X-ray confirmed orbital period, as illus- 
trated by a zoom on two years of data in the left centre panel. Using 
the analysis methods of Section [4] the folded light curve profile is 
parameterised by PS=0.2 and PA=2.5 which is entirely consistent 
with a FRED-like orbital modulation. The light curve for the final 2 
years however shows a much more sinusuoidal profile with a sim- 
ilar amplitude to that seen in the previous years, but with a visibly 
different period and shape (right centre panel of Figure \\3\. While 
it is hard to be certain with just a few cycles of the ~80-90d oscilla- 



tions, the period of these sinusoidal variations may not be constant 
through the two years. 

We interpret these changes in the light curve as a clear indi- 
cation of both orbital-derived periodicities and non-orbital period- 
icities, the latter possibly deriving from aliased stellar pulsations, 
being visible in the same light curve, just as in SXP756. The pres- 
ence of the sinusoidal variations in the last two years of OGLE III 
monitoring make it more difficult to automatically detect the orbital 
variations, but visual inspection of the light curve of the last 2 years 
around each expected periastron passage (shown by the dashed ver- 
tical lines in Figure \Tty finds clear increments in the flux at those 
times, showing that the orbital variations indeed continue. If the 
non-orbital periodicities at ~80-90d period do derive from aliased 
pulsations, we can only state that their true period is below 2 days; 
the degeneracy in the period reconstruction shown in Figure|5]pre- 
vents us from unambiguously estimating the true period. There is 
no indication that the pulsation-derived and orbit-derived variations 
are mutually exclusive, as they co-exist in this light curve, but the 
fact that the two sections of the light curve show such clear changes 
indicates both the transient nature of these phenomena, and the ef- 
fectiveness of light curve shape as a way of distinguishing them. 
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Figure 12. Periodiograms and folded light curve analyses for previously undetected signals in SXP3.34 and SXP893. Based on the folded light curve analysis, 
these two periodicities are unlikely to be due to orbital motion. Dashed horizontal lines on periodograms indicate 99 and 99.9% confidence limits. Dashed 
horizontal lines on folded light curves indicate 10% and 50% of the maximum value, where PS and PA are evaluated (see text for more details). 



5.4 Transient phenomena 

Several of the sources discussed here display modulations that are 
transient in nature. In the case of orbital modulations, this is usually 
modelled by the growth and decline of the decretion disk around the 
Be star: when the disk is large, the orbital modulation is more pro- 
nounced. SXP 18.3 is a good example - no periods may be detected 
when considering the entire light curve, but the orbital modulations 
are clearly detectable in a section encompassing approximately two 
years of monitoring. Other examples may be less extreme but al- 
most all the sources show some variability in the amplitude of the 
orbital modulation. 

The sinusoidal modulations in systems, that we believe to de- 
rive from aliased stellar pulsations, are also often seen to be tran- 
sient in nature. While many are seen at all times throughout the 
OGLE II and III datasets with almost constant amplitude, oth- 
ers are seen to come and go, and even change in frequency (see 
next section). Moreover, we sometimes see evidence of very abrupt 
changes in the pulsation-like behaviour - a classic example being 
given by SXP7.92 in Figure [14] which shows three consecutive 
years of OGLE III data. The year-to-year amplitude shows signif- 
icant changes, but even more dramatic is the abrupt cessation of 
pulsations at TJD~2940. 



It is also appropriate to point out here a limitation of the meth- 
ods used - the Lomb-Scargle method is insensitive to modulations 
superimposed on large non-periodic changes in brightness of the 
underlying star. We have attempted to overcome this as far as pos- 
sible by the use of detrending on appropriate timescales in order 
to restore a level baseline against which the Lomb-Scargle method 
works best. However in some systems, rapid underlying changes 
can dramatically reduce the sensitivity of the period searching for 
some parts of the dataset. Thus it is important to separate out tran- 
sient phenomena intrinsic to the source itself and those introduced 
by changes in sensitivity to modulations. 



5.5 Period drifters 

One feature that has been suggested in past studies of the pulsa- 
tions in Be stars in such X-ray binaries is the phenomenon of pe- 
riod drifting - i.e. slowly varying periods, probably associated with 
Non-Radial Pulsations (NRP) from the Be star l |Coe et al.||2002[ 
|Schmidtke~ & Cowley 20051. Though the expectation for such pe- 
riod changes being associated with NRP behaviour does not com- 
fortably fit with stellar oscillation models, it is extremely unlikely 
that such changes could be associated with a changing orbital pe- 
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riod on the timescales of the OGLE data. |Coe et ah] ( |2002| > dis- 
cussed this possibility in the context of a possible triple system 
for SXP323. Since aliased non-radial pulsations are believed to 
be responsible for the sinusoidal modulations in many OGLE light 
curves, we have investigated this phenomenon of a drifting period 
in several systems where the periodicity is strong enough to deter- 
mine the period accurately on annual timescales. 

• SXP214 shows a clear and steady drift in period - the first 
time a drifting periodicity has been determined in this system (Fig- 
ure [15}. To determine the true drift rate we assume that the mea- 
sured period of around 4.55d shown in the figure is actually the 
beat between the true period and the one day sampling of the OGLE 
data, but as stated previously, there is a degeneracy in this possible 
solution. Assuming a constant drift rate over the period of the ob- 
servations, the degeneracy indicated in Figure|5]leads to two possi- 
ble solutions for the true period and drift: either a period of 0.818d 
with an associated drift of +15.0s/yr or a period of 1.277d with a 
drift of — 37s/yr. 

• SXP323 shows a different behaviour (Figure[T6]l. While there 
is still an overall drift in period, it appears there may not be a sim- 
ple monotonic progression like that seen in SXP214. Nonetheless 
a simple linear fit to the data over the duration of the OGLE ob- 
servations gives a change in the period of +4.92s/yr. Superimposed 
on this steady change is an apparent long term modulation with a 
period ~2000d, but the limited length of the data set means that the 
evidence for such a modulation is tentative at best. 



• SXP701 also exhibits a convincing period drift over the dura- 
tion of the OGLE observations. In this source there is evidence for 
two possible observed periods (-~0.2845d and ~0.282d) with the 
main power jumping from the former to the latter after two years. 
If a linear fit is made to the most persistent second period then a pe- 
riod drift of +5.26s/yr is obtained - similar to that seen in SXP323. 
We note that |Schmidtke & Cowley] ( |2005| > obtained -2.6s/yr for the 
period drift in this object from the shorter MACHO data set - the 
negative rate of change possibly arising from the misleading effect 
seen here in the OGLE data whereby different frequencies can take 
turns to dominate. 

There has been some discussion in the literature as to whether 
the occasionally reported period drifts in Be stars (eg |Stefl &| 
|Balona|(l996| l) are real, or due to unresolved beating between sim- 
ilar periods with differently changing amplitudes. A discussion on 
this topic may be found in Baade (1999) but with no clear con- 
clusions. Certainly the data presented here represent some of the 
longest time bases for such observations and further investigation 
into the OGLE archives for similar stars and periods may be prof- 
itable. 



5.6 Superorbital periods 

It has been suggested by Rajoelimanana, Charles, & Udalski 
< |20 1 1 ) > that these systems may exhibit optical modulation on longer 
timescales than their orbital periods. Furthermore those authors 
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present results which indicate a possible correlation between or- 
bital and superorbital periods (see Figure 45 in their paper). Though 
their plot looks convincing, care must be taken when creating such 
data sets. Since the total length of the OGLE data is ~5000d, then 
any period longer than half of that must be treated with some scep- 
ticism. This effectively removes 6 of the 19 objects from their plot 
at the longer period end. 

At the shorter period end the work presented here casts doubt 
on the value for several periods quoted in Rajoelimana na, Charles,| 
|& Uda lski (20111 as orbital. If the shape of the modulation is 
strongly sinusoidal then there is a very good chance that it arises 
as an alias of an underlying NRP. For example, the period of 36.4d 
given for SXP7.92 appears from our work to be a strong case for 
being the alias of an NRP, as is also the case for SXP101 (21. 9d) 
and SXP3.34 (11. Id) - the latter which Rajoelimanana ,~C"harles, &| 
Udalski (2011) describe as a "marginal" case anyway. If one re- 
moves these objects from their plot, together with the long period 
ones, then the case for a relationship between orbital and possible 
superorbital periods in the remaining 10 objects is greatly reduced. 
There is no doubt that Be stars in the SMC (and probably else- 
where) show considerable creativity in their modes of variability 
(Mennickent et al. 2002), so great care must be taken when claim- 
ing actual periodicities when, instead, such observational changes 
may represent general long-term timescales for change in the mass 
outflow from such stars. If such superorbital periods exist, then the 
warping of the circumstellar disk by an inclined neutron star orbital 
plane may provide an explanation (M artin et al.|201 1) . 

5.7 The V-band light curves 

Although the V-band data coverage is only ~10% of the I-band 
data, we searched the V-band data to determine if any of the 
types of periodicities were more systematically visible in that band. 
A few periodicities were significantly detectable in the V-band 
lightcurves, namely those in SXP7.78, SXP293 and SXP1323; an- 
other, SXP327, is detected but below a formally significant level. 
It is difficult to explain why these particular sources exhibit peri- 
odicities in the V-band while others do not. While we might ex- 
pect persistent pulsations of the donor star to be easier to detect in 
the V-band than any decretion disk modulation, both confirmed or- 
bital and suspected pulsation periods are included in this very small 
subset of detected periods. Although the amplitude of these mod- 
ulations is certainly higher than average, we cannot say for sure 
that large amplitude is the only reason for detection. Unfortunately, 
other sources that we might expect to detect on such a basis, such 
as SXP18.3 and SXP893, exhibit their modulation only transiently, 
and have only 16 and 4 V-band data points respectively in the rel- 
evant time periods. The data for SXP756 show only one periastron 
outburst, which is strong but obviously offers no basis for period 
detection. In all these cases then, we can state that the limited data 
coverage prohibits any possibility to detect the periodicities, and 
we should make no physical inference on the lack of detection in 
the V-band. 



6 CONCLUSIONS 

It is clear that the use of ground-based optical light curves like those 
in the OGLE-II and OGLE-III programs to determine the orbital 
periods of Be/X-ray binary systems can be an effective technique. 
But at the same time, it can be misleading, as other mechanisms 
can produce periodicities that, when combined with typical daily 



sampling, alias into the typical orbital period regime. This can be 
illustrated by use of a P or b-Pspin diagram (Figure \T7) in which we 
plot all the periodic signals detected, regardless of their suspected 
origin. The systems for which the modulation is known, or strongly 
suspected, to be of orbital origin, lie in the traditional zone of the 
diagram associated with the Be/X-ray binaries. However, the sys- 
tems for which the origins of the periodicity may be aliased NRPs 
(or other sinusoidal variations of unknown origin) are outliers to 
this distribution, and coincidentally fall in the region traditionally 
associated with supergiant systems. 

We present a method that allows the shape of the folded light 
curve to be used to distinguish between orbital periods and other 
mechanisms in at least some cases, but we point out that this is 
still not a guarantee. For now, the only unambiguous way of as- 
certaining the orbital periods of the Be/X-ray binaries comes from 
X-ray measurements, either by looking for periodic emission or us- 
ing Doppler methods to fit the orbital parameters directly. 
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Figure 17. A P or b-P S pin diagram constructed assuming all detected periodicities, coded to indicate the suspected mechanisms for the periodicities. 



Schmidtke P. C, Cowley A. P., 2006, AJ, 132, 919 MNRAS, 1085 

Schmidtke P. C, Cowley A. P., 2007, A AS, 38, 724 Trowbridge S., Nowak M. A., Wilms J., 2007, ApJ, 670, 624 

Schmidtke P. C, Cowley A. P., Udalski A., 2006, AJ, 132, 971 

Schmidtke P. C, Cowley A. P., Udalski A., 2009, ATel, 1953, 1 

Schmidtke P. C, Cowley A. P., Udalski A., 2009, ATel, 1992, 1 

Schurch M„ Udalski A., Coe M., 2008, ATel, 1670, 1 

Schurch M. P. E., et al., 2009, MNRAS, 392, 361 

Schurch M. P. E., Coe M. J., McBride V. A., Townsend L. J., Udal- 
ski A., Haberl E, Corbet R. H. D„ 201 1, MNRAS, 412, 391 

Stefl S., Balona L. A., 1996, A&A, 309, 787 

Townsend L. J., Corbet R. H. D., Coe M. J., McBride V. A., Hill 
A. B., Bird A. J., 2009, ATel, 2202, 1 

Townsend L. J., Coe M. J., Corbet R. H. D., Hill A. B„ 2011, 



© 2012 RAS, MNRAS 000,[TlfT6] 



